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The position-momentum quasi-distribution obtained from an Arthurs and Kelly joint measure- 
ment model is used to obtain indirectly an "operational" time-of-arrival (TOA) distribution following 
a quantization procedure proposed by Kochahski and Wodkiewicz [Phys. Rev. A 60, 2689 (1999)]. 
This TOA distribution is not time covariant. The procedure is generalized by using other phase- 
space quasi-distributions, and sufficient conditions are provided for time covariance that limit the 
possible phase-space quasi-distributions essentially to the Wigner function, which, however, provides 
a non-positive TOA quasi-distribution. These problems are remedied with a different quantization 
procedure which, on the other hand, does not guarantee normalization. Finally an Arthurs and 
Kelly measurement model for TOA and energy (valid also for arbitrary conjugate variables when 
one of the variables is bounded from below) is worked out. The marginal TOA distribution so 
^Lj , obtained, a distorted version of Kijowski's distribution, is time covariant, positive, and normalized. 
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00 . !• INTRODUCTION 

00 

In 1965 Arthurs and Kelly [y proposed a model for the simultaneous measurement of the position and momentum 
of a quantum particle. This model has been a valuable tool in many fundamental and applied works, see e.g. |lC| ] 
, and references therein, and is present, albeit implicitly, in many others flT|-|l4||. In the last few years, the advances 
in quantum optics have led to experimental techniques and theoretical approaches devoted to the simultaneous mea- 
surement of conjugate variables that have made possible a practical realization of the original gedanken experiment 
of Arthurs and Kelly §-| , P^|Jl6| . The measurement process is modeled by a Hamiltonian that includes degrees of 
freedom of the particle and two pointers. It is constructed by adding two sudden interaction terms that, when acting 
separately, provide impulsive (or von Neumann) measurements |l7j of the position and the momentum of the particle 
at time t, 

=5 

H AK = 6(t)(7rpp + n x x). (1) 



All other terms are neglected, in particular the ones corresponding to the free dynamics of pointers and particle. In 
each of these two partial interaction Hamiltonians the operator for the particle property to be measured (position x, 
resp. momentum p) is multiplied by the conjugate operator, ttx, resp. 7Tp, of the associated pointer observable fix, 
resp. /ip, so that for each of the ordinary von Neumann measurements the state of the pointer variable is displaced 
proportionally to each eigenvalue Jl7{ . However, the combination of the two interaction terms in a single Hamiltonian 
implies a mutual disturbance of the two measurements. 

In the combined (joint) measurement the displacement of the two commuting pointer variables fix and /ip can be 
described by a true joint distribution, p(/ix, l^p)- It is well known that a unique position-momentum distribution 
for x and p (positive, bilinear in the wave function, and with the correct marginals) cannot be defined in quantum 
mechanics, but there arc many possible quasi-distributions. The Arthurs and Kelly model provides a simple operational 
realization of a family of joint position-momentum quasi-distributions in terms of the commuting pointer positions as 

i 

F(x,p) = p{p-x = x,p, P =p) . (2) 

In fact Arthurs and Kelly considered one particular set of states for the apparatus that makes F a "Husimi function" 
W; however, more general apparatus states are possible. The resulting family of quasi-distributions has been discussed 



by several authors 18 1 



The Arthurs and Kelly process also provides a natural way of quantizing classical functions of position and mo- 
mentum. In particular, for the free motion case, Kochahski and Wodkiewicz have recently defined "operationally" a 
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time-of- arrival (TO A) distribution from the phase-space distribution p |L9]. As in (g), the idea is to use the commuting 
variables px and pp, instead of the non-commuting particle position and momentum, in the classical expression that 
defines the arrival time. 

The theoretical treatment of time observables is an important loose end of the standard quantum mechanical 
formalism. Among these observables, the time of arrival has been investigated using many different approaches, as 
reviewed in - for more recent works see [^T|-^l|. They may be classified according to their intrinsic (ideal) or 
operational nature. Intrinsic models consider the particle on its own, without any external influence other than the 
potential in which it moves, and provide ideal quantities that do not depend on any other degree of freedom. An 
example is the distribution of Kijowski p2| , which satisfies in an unique manner a number of properties motivated by 
its classical analog, or the distribution obtained within the causal theory of Bohm P^| . Operational models, instead, 
include extra degrees of freedom for the measuring device explicitly in the Hamiltonia n [|l9||34|]35[ ] , or implicitly by 
means of a non- unitary evolution law (e.g., with effective complex Hamiltonians |3^j20| |). Notice that the word 
"operational" has a rather different meaning in the usage of some authors J37j], but we shall use it in the sense 
described. This paper investigates the properties (in particular covariance) of operational time-of-arrival distributions 
obtained by means of Arthurs and Kelly measurements, both of position and momentum and of the conjugate variables 
energy and time of arrival. In order to attain a broader perspective, we also consider and put forward other possible 
TOA distributions, based on different phase-space quasi-distributions, not directly derived from an Arthurs and Kelly 
measurement, and look into the relation between operational and ideal quantities. 

In Sec. II we shall for completeness rederive the phase-space distribution p(px , A*p); using the treatment of the 
Arthurs and Kelly model given by Appleby In Sec. Ill we show that the TOA distribution thus obtained, 

following Kochahski and Wodkiewicz, is not covariant. We analyze in the following section alternative distributions, 
with particular emphasis on the property of covariance, and conclude that, within a broad family of phase space 
functions ]3j| , |39| ], only the Wigner function gives a covariant TOA distribution when following the recipe proposed 
in p^. This covariant distribution, however, is not positive, and this defect, in turn, leads us to investigate other 
quantization recipes to obtain covariant and positive distributions, a task in which we succeed, using as before the 
simultaneous Arthurs and Kelly measurement of position and momentum. The result however is not normalizable. 
Yet another option in order to obtain covariant TOA distributions is to perform an Arthurs and Kelly measurement 
for the conjugate variables time of arrival and energy, as we show in Sec. V. This is a simple model for the important 
experimental need to know the energy and the time when the particles arrive at detectors. The results follow closely 
but not exactly the standard Arthurs and Kelly process for position and momentum, and we are again successful in 
obtaining covariant TOA distributions, now correctly normalized. In the final section of conclusions we indicate other 
possible extensions of our results. 



II. THE ARTHURS AND KELLY MEASUREMENT FOR POSITION AND MOMENTUM 



Recently, Appleby has studied thoroughly the concepts of accuracy and disturbance in the Arthurs and Kelly 
model HI?]]. We shall show how to recover, following his analysis, the results for the probability p(px,pp), previously 
obtained in a concise manner by Kochahski and Wodkiewicz. The interest of this detailed rederivation is that the 
meaning of the apparatus dependent "window" function is made explicit. The window function provides /?, when 
convolved with the particle's state, see (|l(J below. Our derivation also applies to a class of apparatus states more 
general than that in p9|. 

The operators for the six variables involved in the Arthurs and Kelly process satisfy the commutation relations 

[x,p\ = \px,n x ] = [pp,n P ] = ih. (3) 

Any other pair of these operators commutes (i.e. pp and p x , x and nx, etc.). The unitary evolution operator 
describing the measurement process is 

U X p = e -*^ p P +7x ^ . (4) 

There is no explicit reference to the measurement instant t here. Uxp may be regarded as the evolution operator 
connecting the states of the system before and after a sudden interaction ([!]). It is also possible to interpret it as the 
evolution operator that gives the final state for an interaction K(npp + ttxx) acting during a time At = 1/K, with 
K sufficiently large so that all other terms in the Hamiltonian can be neglected during the measurement time At [[l] . 

If the initial state of the global particle+apparatus system is given by the product state \ip ® ipap), the probability 
distribution for the result of the measurement takes the form 
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p(p x ,pp) = / dx\(x,p x ,pp\U X p\tp ®1pap)\ 



(5) 

This is the key expression that relates the probability distribution of the pointer variables px and pp to the particle's 
initial state. Let us emphasize that p(px, Pp) is a true probability distribution of its commuting variables, fix and fip; 
however, when interpreted as a function of the quantum mechanical particle variables, namely as F(x,p) = p(px = 
x,pp — p), it is to be understood more properly as a quasi-distribution. This comes about because the marginal 
distributions of F for x and p are not the quantum mechanical ones for the state of the particle, even though they 
are correctly normalized and positive. 

To arrive at an explicit form for such a distribution and to describe the experimental errors in the measurement 
process, Appleby introduces the initial and final "Heisenberg picture" operators Oi — O and Of — U xp OUxp, where 
O can be any of the six operators in (^|). In terms of these operators several "errors" are defined. In particular the 
retrodictive error operators, 



txi — Pxf — Xi , 

epi = ppf — % j (6) 

provide the accuracy with which the result of the measurement reflects the state of the system before the measurement 
was carried out; and the predictive error operators, 



ex/ = pxf - x f , 

e Pf = p Pf - p f , (7) 

give the accuracy with which the result of the measurement reflects the state of the system after the measurement. 
Expressions for the final operators in ^) and (Q), 

Xf = U XP xUxp = X + TTp, 

Pf = U xp pU X P = P-TT X , 

fit - ^ - 

P-Xf = U xp HxUxP = px + X + -7Tp , 

-t - ~ - - 1 - 

HPf = U XP p P Uxp = pp +p - -ttx , (8) 



are easily obtained using the commutation relations (||) and the following relation, valid for two arbitrary operators 
A and B, 



e jAg e - 7 A = g + 7 ^ B] + X- \ A, [A, B] 



7_ 
3! 



-4, 



A, [A,B] 



(9) 



The distribution of measured values (^) can finally be written after some algebra as a convolution in phase space R] 



p(jJ-x,fJ>p) = 



OO rOO 



dxdp W tl (p x - x,p P -p)w(x,p) 



(10) 



— ooJ —oo 



where w(x,p) is the Wigner function of the initial state of the particle (just before the measurement takes place), 



w(x,p) 



dye 



7 I'll 



y 



and the apparatus dependent "window" , or "filter" function in phase space, 



W ei (p X - X,p P -p) 



dye R 



ex. 



y 



tXi 



(11) 



(12) 



where epi = pp — p, and ex% = Px — x, is the Wigner function corresponding to g ei , the reduced initial apparatus 
state density operator for retrodictive errors, 



{e X i\QeA(-Xi 



de X f (e X i, tXf\^a P }(^apW X i> e Xf) 



(13) 



The trace is taken over the predictive error of position. Note that the apparatus state has been represented in the 
basis of the complete set of commuting operators (for the apparatus space) of retrodictive and predictive position 
errors, ex% = Px + t?p/2, and ?xf = fix — np/2 (that refere exclusively to the apparatus, and whose conjugate 
momenta are — epj and ep / respectively) . By chan ging th e distribution of retrodictive errors it is possible to obtain 
a family of operational phase space distributions p,^0|j37| . 
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III. TIME-OF- ARRIVAL DISTRIBUTION: INDIRECT APPROACH 



By assuming a separable, pure state form for g fi — \(f>){(f>\ in (|12|) , and using expression (|ll| ) for w(x,p) in Eq. 
the distribution of measured values of position and momentum can be written as 



p{px,pp) 



dx' <l>*(n x - x')^{x')e- l ^ px ' / h 



(14) 



where the "filter" or "window" function <fi(px — %') is the retrodictive error wave function of the apparatus. (It 
is to be noted that while the form of ( |l4| ) depends on the assumed factorized structure of the reduced retrodictive 
error density operator, (O) is more general and does not require such structure.) Here we may see that the quasi- 
distribution functions F^x^p) = p(px = x,Pp = p) are, in the language of time-frequency analysis [[39]], nothing but 
"spectrograms" , that is, the square modulus of Fourier transforms of the particle state multiplied by an apparatus- 
dependent window function selecting a limited spatial region. Clearly, the transform reflects the combined properties 
of the particle state and the window function. A natural condition to impose on 4>(exi) is that it be centered at 0, so 
that the average retrodictive error of position is zero, and the spatial regions where tp and p have significant values are 
in good agreement, at least on average (since otherwise the pointer would be displaced with respect to the particle's 
position). Similarly, for the momentum representation version of (^) we shall demand that the Fourier transform of 
the window function be also centered at zero retrodictive error of momentum. In summary, and using the language 
of |(| , we shall assume that the measurement is "retrodictively unbiased" . 

Note that in the argument of <j) the sign of x is different from the one in jll| . This can be traced back to different sign 
conventions for the momentum pointer and the commutation relation with its conjugate variable. Both conventions 
lead to equivalent results if appropriate sign changes are taken into account. With the present convention, and for 
unbiased measurements, the pointers are located (on average) at the particle's average position and momentum, and 
thus p(px,pp) tracks directly (of course with the unavoidable distortions inherent to the joint measurement) the 
particle's position and momentum distributions. In fl9|j , instead, the position pointer and the particle's position have 
opposite signs when the filter function is centered at 0, see Eq. (21a) of |19| with go = 0. 

Kochahski and Wodkiewicz have used the operational phase space distribution to define, indirectly, a TOA distri- 
bution for the free particle Jl9] | . The basic idea is to work with the commuting operators fix and pp instead of x 
and p. A classical particle with position x and momentum p at time t takes a time T = —mx/p, measured from t, to 
arrive at the origin. This motivates the definition of a distribution of arrival times by means of the following average 
in the "observed" phase space px,pp, 



n 



K W 



(T;t) 



mpx 
p P 



T 



dpx dpp 5 



mp x 
Pp 



+ T) p(p x ,pp,t) , 



(15) 



where —mpx/pp is used instead of the classical expression. (Note that T is the time interval from t to the arrival 
instant t + T.) 

As an example, and in order to have a better grasp of the properties of this distribution, assume, as in |19l, that 
the particle initial state is a Gaussian wave function with initial mean position xq, initial mean momentum hko and 
position "width" 5 (square root of the variance) , 



ip{x,t) = 



2d- 



1/4 



4fe^ 2 



S 2 



•inn 



(16) 



and that the initial "apparatus filter-function" is given by an unbiased Gaussian with width ct, 

(t>{p x - x) 

For the example at hand the distribution p(px, pp,t) can be written explicitly: 



9 \ V 4 
_f_) e -(p*-*) 2 /- 2 

7TC7 2 / 



(17) 



p(px,pp,t) 



aS 



1 



2nh 



: exp 



1 S 2 a 2 (6 2 + <7 2 ){p P - hk a ) 2 



2h z 



x exp 



^l^f + {5 2 +a 2 )- 

hkpt 



^f + (S 2 + a 2 ) 2 



4t 2 h 2 




(18) 
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Inserting this expression into (|T^) we obtain 



n ^ (T;<) - 2^ MTJ) P { W+(^+ ° 2 ) 2 ) 

l + yfrt(T,t)e?V>H[ttT,t)]) . (19) 



where $ is the error function and the following symbols have been used: 

A(T, t) = ^£±pL + + lM* + ^ (20) 

m 2 m 2 4 

, koth 

Td(t) = - °, m , (21) 



o 



M t) = 2fc [A(0, 0) + H(^(T + TOO) + S 2 TT cl (t))/m] _ 



V2A(T, i)^/(<5 2 + o- 2 ) 2 + 4t 2 h 2 /m 2 



This distribution is shown in Figure 1, where its lack of covariance is evident. According to the expression (|l9|), the 
probability for arriving at instant T + 1, when an interval T has passed after the reference time t, is not equal to the 
probability for arriving at the same instant (T + 1), when an interval T — t' has passed after the reference time t + t' , 

U KW (T;t)^n KW (T-t';t + t'). (23) 

Covariance is however a basic physical requirement for any good quantum time-of-arrival distribution fil] J^^H , It 
simply means that the number of arrivals predicted for a particular fixed instant (T + 1 in p3| ) should be a constant 
quantity independent of t', i.e., on the reference time used in making the prediction. A good apparatus should give a 
stable, fixed answer, independent of the instant that we have switched it on. Lack of covariance implies that different 
predictions are given about the number of arrivals at the same instant of time depending on the reference time chosen 
to make the question (that corresponds here to the Arthurs and Kelly measurement). We may wonder whether a 
time distribution obtained from other phase-space representations of the particle state can be covariant with respect 
to time translations. This study is carried out in the next section. 

IV. COVARIANCE IN TIME 

The question we shall address first is whether any of the many possible representations of the quantum state, other 
than the spectrogram, does provide a covariant time-of-arrival distribution following the recipe shown in ([l5|). The 
process of Arthurs and Kelly associates states of a particle with probability distributions on the phase space of the 
particle, via Eq. (g). This association, in the case of Arthurs and Kelly measurements, is done through the filter 
function, derived from the state of the measuring apparatus. There are however many other ways of building up such 
a pairing between particle state and particle phase space. (Only some of them have a simple operational interpretation 
in terms of a measurement model). 

A very broad class of quantum quasi-probability distributions F of position and momentum was studied and defined 
by Cohen p8|,[39| , the one in Eq. (|l4|) being a particular case (see below) . In Cohen's approach each of the distributions 



F is obtained with a different kernel % from the density opera tor g of the particle, see Eq.(Al) in the appendix; and 



for each of these kernels a quantization rule is defined, Eq. (A4), such that the expectation values can be equally 



obtained by means of phase space integrals or operator traces, Eq. (A3). In particular, 



U 5 (T;t; [ X ]) = (s(j^+t\\= J°° J°° dxdp s{^j+t\ F(x,p,t; [ X ]) - Tr[g(t) S X (T)\ , (24) 

where S X (T) is a shorthand notation for the operator corresponding to the classical "function" 5(mx/p + T) by means 
of the "x-quantization rule" , 



S X (T) = ^I I I I dxdpd0dT8(T + zm/p)x(O,T)exp[-i(O(x-x) + T(jp-p)]. (25) 



OO POD POO pOO 



oo J — oo <J — ocJ — OO 



Covariance in time of the arrival time distribution means that II^ (T; t; [%]) should be equal to 
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n s (T -t';t + t'; [x]) = Tr[g(t + t')S x (T - t') 



Tr 



g{ty St ' /h S x (T - t')e- tSt 'l n 



for all g(t). Thus the covariance condition can be expressed in operator form, 

S X (T) = e ig *'/ a ? x (T - t')e~ lSt ' /h . 

Our next task is to find kernels x that fulfill (^8|) . To this end we shall work out the momentum representation of the 
two sides of (Eq) for an arbitrary \. For the left hand side we find: 



(26) 
(27) 

(28) 



(p'\S x (T)\p") = 



1 



oo poo 



dpdr\p\x 



oo J — oo 



p -p 



2-Khm 

e i(p' 2 -p"' 2 )T/(2mh) /-oo j-co 



e ~ip(r-T(p ~p )/mh) e ir{p+p )/2 



27rhm 

Use has been made of the relation 



dpda\p\ x 



oo J — OO 



p -p 



T(p'-p") 



>-[p+(p'+p")/2] 



(29) 
(30) 



e i(9q+rp) _ e ih9r/2 e i9^ e . 



and of the change of variable a = 
hand side of (G8f) takes the form 



a*TP (31) 

r — T(p' — p")/mh. Operating similarly, the momentum representation of the right 



(p'\e lSt '/%(T - t')e- iSt '/ h \p") 



3 i(p' 2 -p" 2 )T/(2mh) 

2Tthm 



dpda\p\ x 



OO J — OO 



p - p 



(T-t'W-p") 



Jcr[p+(p'+p")/2] 



(32) 



If ( p2[ ) has to be equal to ( p9| ) for all t' , the kernel function x(0,r) must be independent of its second argument. 
Since the only set of kernels worth considering are those that preserve the normalization of the state, i.e., such that 
x(0, 0) = 1 1 38 3S|] , independence of r implies x(0? t) = 1, which is a sufficient condition to provide a correct marginal 
distribution for p. Therefore, independence on r limits the set of possible x kernels rather strongly, in particular the 
spectrogram does not generally belong to this class, since it does not generally satisfy the marginals. Its kernel has 
the following form, see (14), 



X(fi,T) 



dy<j)*(y - Th/2)(f>(y + Th/2)e iev . 



(33) 



This may become independent of r in the limit of a very flat <f> function, namely, for a vanishing retrodictive error 
of momentum. This is not a desirable limit for our purposes: the position becomes so imprecise, that Hkw(T) also 
tends to vanish in that limit, as can be seen, for example, by taking a — * oo in Eq. (|l9|). 

Restrictingourselves to the x functions that provide an F with the two correct marginals (this requires x(@> 0) = 
x(0, r) = 1 |3^j39[| , and excludes the spectrograms), the independence on r can only be satisfied by x — 1> which is 
the kernel that corresponds to the Wigner function and the associated Weyl quantization rule. Furthermore, in the 
set of the kernels which provide scale-invariance ]39| l, characterized by being functions of the product t0, x = 1 is 
again the only possible case that presents covariance. 

The covariant TOA distribution obtained with the Wigner function and (p3) takes the form 



1 />oo /-OO 
Rs,Wigner{T;t) = — / / dp' dp 

tim J _ oc J _ oc 



P +P 



(p'\g(t + T)\p") 



-ooJ — oo 

= Hs,Wtgner{T — t';t + t') . (34) 

Let us compare this result with the flux at the origin, J — dpw(0,p, By substituting (|ll|), and performing 
several integrals, this takes the form 



J = 



1 

hm 



oo c oo 



dp'dp" 



OO J — OO 



p' + p h 



(p'ldW 



(35) 



The only difference with (|3J) is the presence or absence of the absolute value in the half-sum, the two results being 
equal for states without negative momentum components. However, a quantum mechanical state composed by positive 
momenta is compatible with a negative value of J (backflow) at certain times and positions |^,[l3|,^7]] . Therefore, 
the "time-of-arrival quasi-distribution" obtained by means of Wigner's function satisfies the covariance under time 
translations, but fails to be a positive distribution, which is also an important requirement for a good time-of-arrival 
distribution. 



G 



A. Kijowski's distribution 



It is also of interest to investigate whether there is a kernel function x such that the corresponding distribution 
IT(T) obtained via is equal to the covariant distribution of Kijowski, 



U K (T;t)=Ti 



g(t) i^2\T,a){T, 



Here \T,a) are the (generalized) eigenfunctions of the Aharonov-Bohm time-of-arrival operator (see below), 

\P\ ^ V2 



(p\T,a) 



Jp*T/2mha 



(36) 



(37) 



and a — ± is the degeneracy index associated with positive or negative momentum. 

We look for a kernel x such that S X (T) — ^2 a |T, a){T, a\. Working again in momentum representation, 



(p'\8 x (T)\p") = Q{p>p") 



V'W 2 



(P'P") 



mh 



e i(p -V ,2 )T/(2 m h) e i(p'-p")X/h 



(38) 



Comparing with (|29|) we sec that x need not depend on its second argument, as should have been expected since Hk 
is covariant. The integrals in (|2^) can be then carried out, and the requirement (|38l) becomes 



X 



h 



&(p'p"Wp") 1/2 



p 1 + p" 



(39) 



Changing variables to v = p' — p" and r\ = (p' +p")/2, it is clear that no function of v can satisfy this equation, since 
the right hand side also depends on n. It follows that Kijowski's distribution cannot be obtained by such an extension 
of the procedure of Kochahski and Wodkiewicz. 



B. Alternative quantization of time of arrival 



Up to now we have used position-momentum representations of the particle state to define time distributions by 
means of (^4|), i.e., by an extension of the procedure proposed by Kochahski and Wodkiewicz. Let us now show that 
this does not exhaust all possible quantizations of the classical distribution, and in fact that it is possible to do much 
better with respect to covariance and positivity. Consider the following quantities, 



nj(T;t; [ X ]) 



dx 



dp 



5(x)F(x,p;T + t; [ X ]) , 



(40) 



where the subscript J indicates that in the classical version of this expression, where F would be a true phase space joint 
distribution, ILj represents the positive minus the negative fluxes (the positive flux is defined by J dp F(0,p)p/m, 

and the negative flux by J_ dpF(0,p)p/m). This is apparently quite different from (pi| ) but, in fact, in the context 
of classical mechanics of the free particle, the two expressions are equivalent, as can be seen by using Liouville's 
theorem and the trajectory equation for free motion. However, their quantizations are not the same in general. Note 
in particular that in (p^j) the time dependence has been put entirely in the state and not in the observable. An 



important consequence is that the functions Hj(T;t; [x]) are time covariant for all x- A case where ( p!4[ ) and ( |40) ) 
are equal is X— 1 ; because the propagator of the Wigner function for free motion is just the classical propagator. 
In general, (ff0j) can be understood as a "covariantization" of (U): they coincide for T = 0, and (f4(i| ) is covariant by 
construction, which means that 



Uj(T;t; [ X }) = n s (0:t + T; [ X ]) 



(41) 



for all t and T, and function x- 

The difference between ( |40| ) and ( p4[ ) is actually more profound: whereas in ( p4| ) we are performing quantization of 
a family of classical functions, parameterised by T, namely, 6(mx/p + T), ( |40| ) pertains to a family of quantizations 
for the same classical function, |™|^(^)) each of the quantizations being parameterised by T. Since the value of 
Kijowski's distribution for a given T is actually the trace with the state of a suitable operator, it seems unlikely that 
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it can be reproduced within a scheme so radically different in conceptual terms as a family of quantizations, and, 
in fact, using a procedure similar to the one in the previous subsection, one can prove that Kijowski's distribution 
cannot be obtained from the quantization of |^| S(q) with any x within the class being considered 0]. 

Even so, positive TOA distributions can be obtained from expression ( f!o|) for suitable quantization functions x- 
Indeed, a sufficient condition is that F[x] be positive. In particular, the spectrogram is positive and provides a 
covariant TOA distribution via (fiol). For the example worked out before, see ([l8|), 



h Sa^-^P^ + (S^ + a^ ( -2k^T 2 (t + T) + a^(0) + A(0,0)A 
U ^ i,t; -27rm A(0,i + T) ° Xp ^ + ((?2 + a2)2 j 

l + V<(0;T + i)e^ 0;T+t )$[£(0;T + t)]) . (42) 

This distribution is represented in Figure 1. The problem is that, due to the asymptotic dependence ~ 1/|T| for large 
\T\, it cannot be normalized. In fact whereas by construction (|24|) is normalized (provided the phase space density is 
normalized), (|4^) is not automatically normalized for an arbitrary Xi an d can actually be non-normalizable as in this 
example. Positive, covariant, and normalized operational TOA distributions will be obtained in the following section. 



V. ARTHURS AND KELLY MODEL FOR TIME OF ARRIVAL AND ENERGY 

In the previous section TOA distributions have been defined indirectly from position-momentum quasi-distributions. 
Similar energy distributions can also be obtained indirectly, as pointed out in p9| . However, an Arthurs and Kelly 
type of measurement can be used to find operationally TOA and energy distributions in a direct way: instead of 
the conjugate variables position and momentum, we will describe the measurement process in terms of the conjugate 
variables energy E and time of arrival T. (In this section we shall use the somewhat unusual notation E for the free 
motion Hamiltonian of the particle E = H = p 2 /(2m) to keep the notation in parallel with the {x,p} measurement 
model.) In our {E,T} model the operator T is taken as the time operator introduced by Aharonov and Bohm [Q, 

to Al 1 A 
2 V P P ) 

by symmetrizing the classical expression, —mx/p, for the time of arrival at X = 0, computed from t, of a particle that 
at time t has position x and momentum p. This operator is not self-adjoint but maximally symmetric, see a detailed 
discussion in p8[ . Its (generalized) eigenfunctions in momentum representation are given by (|37|). The {E, T} Arthurs 
and Kelly measurement process is based on pointer and particle operators, JLtiV'Ei^Ti^Ei E, and T, parallel to the 
set used for x and p, and related by similar commutation relations, 

[E, f] = [p E ,K E ] = \MT,*T] = ih, (44) 

with all other commutators being zero. 

The energy operator is bounded from below. As explained later it is also necessary to require that the operator 
7Tt be bounded from below. Another difference with the {x,p} case is the degeneracy, associated with positive and 
negative momenta, of the spectra of E and T. We shall keep the pointer variables non-degenerate for simplicity 
although a more detailed model including pointers sensitive to the degeneracy index is also possible. 

The evolution operator describing the measurement process is now given by 

U ET = e-^T+^£) _ ( 45 ) 
Our purpose is yet again to obtain the probability distribution for the result of the measurement, 

/•OO 

P{He,V>t) = / dE^E.a^E^rWET^^^ap)] 2 ■ (46) 

ol J 

In Sec. II, in the {x,p} phase space, the error operators and all their relations were obtained using the definitions 
of final and initial operators, the commutation relations (|^), and expression d). In the {E, T} phase space similar 
relations are also valid, and (0) to rt|) hold by substituting E for x and T for p. 



8 



In order to construct a meaningful Arthurs and Kelly model for TOA and energy, it is useful to analyze first the 
simple, but not trivial, von Neumann measurement model of T, corresponding to the operator exp — vkiT/U. Note 
that the basis {\T, a)} is complete but non-orthogonal, because T is not self-adjoint. However, the intermediate com- 
putations are mostly carried out in the H-E^a)} basis, which is indeed complete and orthogonal, since it corresponds 
to the spectral decomposition of a self-adjoint operator. Using the overlap 

(E,a\T,af) = -^e iET ' n 5 act , , (47) 

we find 

(E, a, TT T \e~^ TT E', a', tt' t ) = 5{E - tt t - E')5 aa ,5(-K T ~ n' T ) . (48) 

When integrating over E, which is a positive variable, the energy delta function cannot be satisfied unless E' + ttt is 
positive. This means that the exponential displaces the energy of the energy eigenstate as long as E' + ttt remains 
positive, but annihilates the state otherwise, 

e -^ TT \E',a') = G(E' + Tt T )\E' +ir T ,a') . (49) 

It would be more rigurous to deal with wave packets in the energy representation, given that T is not self-adjoint. 
However, the results we will be using are actually unchanged, so there is no need to insist on this fact. In order to 
guarantee that the exponential operator is isometric (i.e. preserves norm), we shall restrict the spectrum of ttt to the 
positive half line, ttt and ht are thus a conjugate pair of operators similar to E and T except for the degeneracy. 



In particular, this means that ht is not self-adjoint, because of Pauli's theorem [J46 28|. Using similar arguments for 
each of the exponentials in 

Uet = e -i(^ T +^ s ) = g-^g-^e^^ , (50) 

it is easily proved that Uet is isometric provided that ttt is positive, and irrespective of the (real) spectrum chosen for 
He- Indeed, there would be two different possible models, depending on whether < he < oo or — oo < He < oo is 
chosen. The bounded case has the advantage of strictly positive values, mimicking the restriction on values of E, but 
the disadvantage that initial apparatus states near the origin cannot be Gaussian nor symmetrically peaked around 
He = 0. For the unbounded case, the initial state may be located symmetrically around the origin, however, the 
price of there being negative values has to be paid. We shall work out the second possibility hereafter, but, with little 
changes, the analysis of the alternative apparatus may be also carried out. In particular, the important property of 
time covariance is valid for both cases. Let us first see how Uet acts on a convenient basis for the composite system: 

ttt \ 

E + tt t , a, he + E + — ,ir T ) ■ (51) 



U E t\E,oi,ij,e,itt) 
By inserting a resolution of the identity in this basis one finds 
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(E,a,(j,E,TTT\UET\i> ® Vv) = ( E - n T ,a,HE - E + "ipT 



il)®ip ap j®{E-Tz T ), (52) 



whose Fourier transform (to change from ttt to the conjugate variable he) is, using now the basis of retrodictive and 
predictive energy errors for the apparatus state, and E' = E — ttt, 

(E 7 a,HE,HT\UET\*P®*Pa P ) = J^ I dE'e^^-^ifiE - E',he - E\yj ap )(E',a\^} . (53) 

The upper limit reflects the fact that the measurement always shifts the particle's energy "upwards" . Any resulting 
energy component (E\ can only have contributions from the lower energy states. From here the operational joint 
time-of-arrival and energy distribution takes the form 

-i />oo r-E r-E 

p[ t i E ,fiT;Ht = 0)} = TJ2 dE dE ' dE"e^E"-E') 
h a Jo Jo Jo 

x {he - E',he - E\^ ap )(^ ap \HE ~ E", he - E)(E', a\^(0))(^(0)\E" , a) . (54) 
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Let us now examine some properties of this distribution: It is time covariant as can be easily seen from (|54|): 

i pOO pE pE 

P \pL E , l x T -f^{t)]=Y J T dE dE' d E"e^-^ E "- E '^ 

a 11 JO JO Jo 

x (fi B - E', (IE - E\^ ap ) {i> ap \ m - E", fi E - E) (E 1 , a\e- iE,t ^(0)) {^(Q)\e iE " ^E" , a) 

= p[M£,Mt;-0(O)] , (55) 

because of the cancellation of the t-dependent exponentials. It is also possible to write it in a form similar, but not 
identical, to Eq. (|l0|). To this end we shall define the following apparatus dependent object in the basis of retrodictive 
error of energy, 



(a\g ap {c)\b) = / dE(a,c-E\ip ap )(ip ap \b,c-E)Q(E-c+a)G(E-c + b). (56) 
Jo 

Changing variables to the half-sum and difference, 

E' + E" 

V = E" - E' , (57) 

we can write (|54|) as 

/>OC poo 

p(he,Ht) = Y i -t ds dye^y^inE-s + y^pinE^E-s-y^) 
a n Jo J-oc 

x(s- y/2, a\fl>) (#> + y/2, a) G(s - y/2)Q(s + y/2) . (58) 



Alternatively, 



pOO poo 

p{he,Pt) =^2 dE dTWfj, E (fi E - E,p, T -T)w a . a (E,T) , 

„ Jo J-oo 



(59) 



where 



and 



1 f°° 

w a , a ,(E,T) = ~ / dyj^iE - y/2,a\i>){i>\E + y/2,a')e(E - y/2)Q(E + y/2) , (60) 



1 f°° 

W^ E (p E -E,p T -T) = - J dy e -<^~ T ^/ h (p E - E - y/2\g ap (p E )\p E - E + y/2) , (61) 

as can be checked by substitution of (|60| ) and (^l|) into ([59|), and renaming s = E in (|58|). is an apparatus 

dependent window function for the energy-TOA distribution of the particle. Even though ( |6l| ) and ( [56] ) provide an 
explicit expression, its interpretation is not as simple as the corresponding phase space window function for the {x,p} 
case of Sec. II. Because of the lower energy bound of the energy, g ap (p E ) is not, in general, the reduced density 
operator of the apparatus for the retrodictive errors. The filtered function, w a . a (E, T), has a simpler content as an 
a, a diagonal component of the energy-TOA Wigner matrix for the particle, Eq. (|60|). Remarkably, no interference 
term with a ^ a' contributes to ([59]), a feature shared with Kijowski's time-of-arrival distribution. As a matter of 
fact, by tracing over a and integrating over E, the marginal of the particle's Wigner matrix is nothing but Kijowski's 
distribution, 

/•OO -| poo />oo 

V/ dEw aa {E,T) = jY^ dE 1 dE"^ E "- E 'y h {E\a\iP){iP\E",a) 
a Jo n a Jo Jo 

= Y J {T,a\^){^\T,a) =Il K (T), (62) 

a 

as is readily seen by undoing the change of variables displayed in ( |57j ) and using Eq. (p!7j). However, whereas in 
the von Neumann measurement the resulting distribution P(ht) is a smoothed version of Kijowski's distribution, 
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see Appendix B, the TO A marginal of the Arthurs and Kelly distribution obtained by integrating ( fffi ) over not 
only smooths but also distorts Kijowski's distribution. This is due to the double dependence on \xe of the apparatus 
dependent window function W^ E (/iE — E, //y — T): aside from w a . a (E,T) the integral over he leaves an extra in- 
dependent function. This is one further peculiarity of the {E,T} Arthurs and Kelly model with respect to the {x,p} 
model of Sec. II. In the later, there is no double dependence on \xx in the apparatus window Wigner function, and 
the marginal is simply a smoothed version of the quantum mechanical distribution of the particle's position, without 
additional distortion. 



VI. DISCUSSION 

In this work we have explored different ways to obtain time-of-arrival (TOA) distributions "operationally", i.e., by 
means of models that include the particle and additional degrees of freedom acting as meters or sensors of the particle. 
An indirect route is to define the TOA distribution from a distribution of position and momentum that is determined 
operationally. The Arthurs and Kelly measurement model is possibly the simplest model where such distributions 
may be obtained. We have shown however that in order to obtain a time covariant result the step from the phase space 
distribution to the TOA distribution is delicate, even crucial. Furthermore, when we consider generic distributions 
of position and momentum the essential property of time covariance proves to be rather elusive. In particular, we 
have proved that for the choice of classical phase space function S(T + mq/p) (that is, the choice of Kochahski and 
Wodkiewicz, [fT9|) there is no distribution of position and momentum within a wide class for which the associated 
TOA distribution is covariant in time. Moreover, we have given yet another instance of the well known fact that two 
classically equivalent expressions may give different quantum distributions, which, in the case at hand, are covariant 
or not. Thus, we have proposed a different classical phase space function, 6(q)\p/m\, for which covariance is always 
granted for a wide class of phase space families of distributions parameterized by T, compare ( 
are the quantum distributions numerically different from the classical ones, conceptually they are also worlds apart: 
whereas the quantization of Kochahski and Wodkiewicz should be more properly understood as the quantization, 
within the same quantization scheme, of a family of classical functions, our proposal is rather a family of quantization 
schemes of the same classical function. 

The question now to be posed is whether any of these different indirect methods of obtaining TOA distributions can 
be selected as being "better" in some suitable sense. From the theoretical point of view, covariance and positivity seem 
to be minimal requirements, therefore selecting positive flux {5{q)\p/m\) over trajectory identification (<5(T + mq/p)). 

Nonetheless, in the case of the free particle, not even the quantization of positive flux provides us with the "best" 
TOA distribution, in that it may be not normalizable, as we have shown with an example, and that for any quantization 
of positive flux the variance of the distribution will be larger than for Kijowski's distribution, which we have shown 
does not lie within a wide class of quantizations of 5{q)\p/m\. We are in this manner led to propose a different 
operational approach, in which time of arrival and energy are directly the "measured" variables. In this case time 
covariance is automatically satisfied. Not only that: Kijowski's TOA distribution can now be understood as a 
marginal distribution deriving from the Wigner function, suitably generalized for the pair of variables energy and 
time of arrival. Operationally, a simple TOA von Neumann measurement provides a smoothed Kijowski's distribution 
whereas the marginal of the time-energy distribution that results from the Arthurs and Kelly process is distorted 
beyond a simple smoothing. In both cases and for a generic apparatus state we see that the variance of the TOA 
distributions hence derived is bigger than for Kijowski's, as was only to be expected from the axiomatic derivation of 
Kijowski's distribution, which selects Hk as the one with smallest variance |32j. 

Additional to our central objective, namely, the study of operational definitions of TOA distributions and their 
relation to ideal measurements of time, we have also obtained a general result concerning operational definitions of 
measurements of pairs of conjugate variables when one of them has a bounded spectrum from below. For instance, 
x and p in the half line. The naive construction of Wigner's function for such a case is actually valid, even though 
p is not self-adjoint. The reason for the validity of the naive expression is, however, far from simple. Since p is not 
self-adjoint and admits no self-adjoint extension on the half-line (technically, it is a maximally symmetric operator 
with (1,0) defect indices), there is no spectral decomposition available for it. Even so, we do have the next best 
thing, namely, a positive operator valued measure (POVM) J37],^8|], and this allows the construction of the basis of 
operators required for the definition of Wigner's function. Moreover, the probability density obtained as the marginal 
distribution for the naively defined Wigner function can be used to reconstruct the whole POVM. To the best of 
our knowledge, the applicability of the Wigner function formalism and its generalizations has not been discussed 
heretofore in the literature for this case. 

On a more general note, we would like to point out that even though the world of "operational models" is in principle 
between the ideal theories (depending only on the particle's state) and actual experiments, the connections are not 



24) and 
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always explicit or obvious. We have in particular emphasized that the phase space probability obtained operationally 
by means of the position-momentum Arthurs and Kelly process may be regarded as one particular (ideal) joint quasi- 
distribution. Many others exist, in particular the ones in the family defined by Cohen, but their links with a particular 
operational procedure are not generally as direct. The fact that Kijowski's ideal TOA distribution, the closest object 
to the classical TOA distribution considering the common properties satisfied, is not obtained with any of these phase 
space quasi-distributions indicates that there is a limit to the flexibility of the class of quantization rules included in 
Cohen's formalism, a warning flag worth considering when dealing with other quantization problems. 

The other connection, between operational results an actual experiments is sometimes taken for granted, but in 
fact it is far from being easily realized. The conjugate variables measured by "Arthurs and Kelly experiments" in 
quantum optics are the two quadrature components of the electric field strength of a single-mode radiation field, not 
the position and momentum of a particle. Nevertheless we expect a similar experiment with particles to be feasible 
soon (see ^7|). In future work we also intend to relate the ideal TOA distributions with operational models that 
describe a continuous (rather than instantaneous) measurement process, thus being closer to an actual time-of-flight 
experiment as emphasized by Aharonov et al. |p4| . 

We thank Rick Leavens and D. Alonso for useful discussions, and acknowledge support by Gobicrno Autonomo de 
Canarias (PB/95), MEC (PB97-1482), CERION, and The University of the Basque Country (project UPV 063.310- 
EB187/98) 

APPENDIX A: BASIC COHEN'S PHASE SPACE FORMALISM 

Given an operator G and a density operator g their phase space representatives 

Q 



-y poo poo poo i y^. 

^ />oc poc poo I ^ 

9M= ^lj 6 Lj t Lj u \ u --2 



G 



_^ e -;[<?(*-«)+^>] x (0 )T ) (ai) 

h \ i[e( q -u)+T P ] 

u + — Tir^— A2 

2 / x{Q,r) 



(G(x,p))=tv(gG)= dx d P F(x,p)g(x,p). (A3) 



are chosen such that 



Conversely, given two space space representatives for the state and dynamical variable, their operators are obtained 
by 

+■ pOO p OC pOO pOO ^ 

Q=7r- \ dx dp d6 dTF{x,p) X (0,r)- 1 e l ^ x - x '> +T(p - p ^ (A4) 

J — OC J— oo <i — oc J— oo 

-i pOO pOC pOO pOO ^ 

G(x,p) = — dx dp d9 dTg{x 1 p)e- l{ex+T ^x{9 1 T)e^ 6x+Tp \ (A5) 

J —OO J — OO J —GO J — OO 



APPENDIX B: VON NEUMANN MEASUREMENT OF T 

Using ( f49| ) the amplitude that results from the von Neumann measurement of the time operator T is given, for a 
semibounded ttt, by 

(T,a, MT |e-" T?/ V(u®?M = (»T-T\^ ap ){T,a\i>(t)). (Bl) 

Since there is only one pointer the wave function of the apparatus, ip a p, depends only on one variable here. The 
probability to find fir in a measurement carried out at t is a smoothed version of Kijowski's distribution, see ([56]), 

/oo 
dT\(n T -T\i> ap )\ 2 II K (T;t), (B2) 
-oo 

which is covariant [P(ht — t';t + £') = P(ht', t)] because of the covariance of IIk. 
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FIGURE CAPTIONS: 

Figure 1: Ukw{T — t; t) versus T: t = —0.2 (long dashed line); t = —0.1 (short dashed line). IIj(T — t; t) versus 
T for any t, (solid line). In all cases x = —2.5, k = 10, a = S = 0.1, and h = m = 1 (atomic units). The bumps on 
the left/right are essentially due to contributions from negative/positive momenta. 
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